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PRECkDli'O  FA3S 


INTRODUCTION 


Dynamical  perturbations  o£  the  auroral  zone  atmosphere  by  au¬ 
roral  activity  have  been  observed  at  large  distances  from  the  auroral 
zone  as  low  frequency  traveling  ionospheric  disturbances  (TIDs), 

(Davis,  1971)  and  locally  as  auroral  infrasonic  signals  (AIWs)  (Wilson, 

1975).  Theoretical  descriptions  of  the  forementioned  effects  have 
been  considered  by  Chimonas  and  Hines  (1970),  Chimonas  (1970)  and 
Chimonas  and  Peltier  (1970).  Additional  analyses  emphasizing  various 
aspects  of  the  atmosphere's  dynamical  response  to  auroral  activity 
have  been  carried  out  by  Blumen  and  Hendl  (1969),  Testud  (1970), 

Francis  (1974),  Chiu  (1976),  and  Richmond  and  Matsushita  (1975).  Most 
of  these  authors  have  concentrated  on  the  far-field  disturbance  with 
the  exception  of  Chimonas  (1970),  who  dealt  with  the  infrasonic  wave 
field  near  a  pulsating  aurora,  and  Richmond  and  Matsushita  (1975)  who 
were  primarily  concerned  with  modeling  the  global  scale  thermospheric 
(  >80  km  altitude)  response  to  magnetic  substorms.  This  paper  describes 
a  theoretical  investigation  of  atmospheric  dynamics  in  the  immediate  vicinity 
of  a  transient  auroral  electrojet.  The  method  of  analysis  is  similar 
to  that  used  by  Chimonas  and  Hines  (1970)  in  a  study  related  to  the 
interpretation  of  midlatitude  TIDs  except  that  the  far-field  approximation 
is  not  invoked.  The  results,  obtained  by  a  combination  of  analytical 
and  numerical  methods,  are  used  to  estimate  the  magnitude  of  the 
pressure,  velocity  and  velocity  shear  perturbations  in  the  atmosphere 
below  the  electro  jet. 


THEORY 

Following  Chimonas  and  Hines  (1970),  the  atmosphere  is  treated 
as  a  compressible  fluid  upon  which  the  electrojet  acts  as  a  perturbing 
force.  Although  both  Joule  heating  and  Lorentz  forces  contribute  to 
the  latter,  it  will  be  assumed,  on  the  basis  of  several  recent  studies 
(Hunsucker,  1977,  Brekke,  1978),  that  the  Lorentz  source  is  dominant. 
The  atmosphere  is  assumed  isothermal,  windless  and  nonviscous. 
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The  standard  fluid  equations  that  are  invoked  for  the  present 
analysis  are  continuity 


.  +  7  •  (PV)  =  0  (1) 

at 


Momentum 


P  +  PV  ♦  7  V  -  p  g  +  7p  =  pF  (2) 

at 


and  heating  (neglecting  the  Joule  source) 


iL- (pp“Y)  +  y.  7  (pp’Y)  =  0  (3) 

at 

where  is  the  mass  density,  p  is  the  pressure,  V  is  the  fluid 
velocity,  g  is  the  gravitational  acceleration,  and  y  is  the  ratio 
of  specific  heats.  The  above  notation  has  been  chosen  to  conform  to 
that  used  by  Chimonas  and  Hines  (1970)  .  The  Lorentz  force  per 
unit  mass  F*  JxB  /fl,  where  J  is  the  electrojet  current  and  B 
is  the  magnetic  field,  is  communicated  by  the  moving  ions  to  the  neu¬ 
tral  fluid  by  collisions.  A  rectangular  coordinate  system  is  chosen 
in  which  the  z  axis  is  vertical,  x  is  in  the  north-south  direction, 
and  y  is  in  the  east- west  direction. 


Under  the  assumption  that  the  unperturbed  atmospheric  density 

-z/H 

po(z)  and  pressure  pQ  (z)  are  proportional  to  e  -where  K  is  the 
constant  scale  height,  the  fluid  equations  (1)  -  (3)  can  be  linearized 
and  Fourier  transformed  in  time  to  yield  the  set 

iu>?-  V  /H  +  7  -V  =  0  (4) 

z 

.  '  i®  V  +  g  H  7px  +  (pL  -px)  g  =  F  (5) 

iuip.  -  V  /H  +  Y-.7  •  V  *  0  (6) 

X  z 

where  the  —  over  the  symbol  denotes  the  transformed  quantity, 

u>  e  iwt  f  (m)  (7) 

and  the  subscripts  1  signify  normalized  perturbations  of  the  form 

<P  ~  P0>. 
po 

Pressure  Perturbation: 

Equations  (4)  -  (6)  combine  to  given  an  equation  for  p^  analogous 
to  eq.  (12)  of  Chimonas  and  Hines: 

f*a  +  “»  _£!l_  -  JL. _  ill 

(u» 2  -  a  *  )  S  (  w  2  -  t»  2)  H  dz 
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(  U)2-  u?g)  4u>*  H2 


where  w  .  =  1/2  is  the  acoustic  cut-off  frequency  and  ua 


v(  y  -  1)  g/  is  the  Brunt  -  Vaisala  frequency.  As  pointed  out  by 
these  authors,  the  substitution 


-  z/2H  t  f 

?1  =  e  4>  2) 


reduces  eq.  (8)  to  the  form 


§*  +  a  +  A ,  4>  *  M, 

a?  1  2  1 


where 


(aiZ-w 


2  2  2 
U>  (»  -#A) 

C2  ( u> 2  -  u»  2  ) 

© 


Y  gH 


Mj  *  e 


-z/2H 
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The  procedure  for  solving  eq.  (10)  used  here  differs  somewhat 
from  that  employed  by  Chimonas  and  Hines.  First,  it  is  assumed 
that  the  electrojet  flows  in  the  east-west  direction  and  that  it  can  be 
described  by  the  current 

J  =  y  Af  (z)  g  (x)  T  (t)  (15) 


In  this  case  the  Lorentz  force  for  a  vertical  B  field  gives  a  form 
for  given  by 


Mx  =  Qj{ z)  G  (x)  T  (oj) 


(16) 


where 

BA 

’  1h 


(17) 


j '  i,)  .  e-j/2H  _ii!L  <18) 

P(z) 


G  (x) 


dg  (x) 
dx 


(19) 


Two  successive  Fourier  transformations  of  eq.  (10)  in  the  spatial 
variables  x  and  z  then  lead  to  the  equation 


(’kx  ‘  A!  kz  +  A2)  ^  =  QwT(^z)  G(fcx)  T  (  m)  (20) 


ii 


Equation  (20)  can*  be  solved  for' an -Infinite  a  tine  sphere,  if  the  functions 
us e'd  in  the  description  of  the  electrojet  current  are  carefully  chosen.  The 
model  selected  for  the  present  Analysis  £.3  shown  in  Fig;  1.  If  a  box  function 
is  selected  to  represent  the  sharp  north- south  boundary  of  an  auroral 


'electrojet: 


g  (*)  * 


'  x‘  <  rj 

fj  <  x  <  ij  +  c 

x  >  a  +  n 


then 


&(x)  * 


dg  (x) 


a  6(x-1])-SCx-7]-9) 


and  the  x  dependence  of  the  right  hand  side  of  eq.  (20)  is  contained 
in  the  transformed  function 


3  =  -i  (e^x  *  -e^x  ( ^  +  Q})  /k:a 


It  follows  that  the  Fourier  inversion  in  x  gives 


0  *  Q<£(k-Z)  T  ( «a)  — L.  £  sin  ui^  (x  -  rj)  -  sin  ujq  (x  -  77  -  c)J 


where 


®o  3  [*A1  +  *2  ]1/2 


If  the  altitude  dependence  of  M.  is  approximated  by 


(*)  *  e 


»z/2K  f ( z) 


1  !\  (0) 


.  z  <  4  (26) 

f  <  z  <  4  +>  \ 

z>  4  +  x 


the  inversion  in  z,  performed  with  the  aid  of  the  convolution  theorem 


(eg-  see  Butkov,  1968)  gives 


-Q  T(«») 

Po  «>) 


/|  +  X 

is  [J 

' 


(Xj)  -  Jo  (X2)]  (27) 


in  which  the  arguments  of  the  zero  order  Bessel  functions  JQ  are 


xi  * 


\j—^~  i_Ai (x  ■  i)2  ■ <*  ■ 8)2 ) 


2  \  1/2 


X2  * 


("Aj  (X  -  ))-  9)2  -  (2  -  »)2) 


A  final  inverse  transformation  is  needed  to  obtain  the  fractional 
pressure  pulse  from  ^  : 


V\  (*»*» 


t)  *  — i—  f  p  eiu,t  d»  * 

1.  1  v^r  /. 


The  inverse  transform  of  the  Bessel  function  in  the  integrand  of  eq. 
(27)  can  be  found  in  tables  for  the  low  frequency  (  ui «  uj^)  limit  since 
the  approximations 


A^  **  - 


-Ai 

c2 
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can  be  made.  The  arguments  (eqs  (28)  and  (29))  then  take  the  form 
b  (ui2-  a2)*^2,  in  which  case  the  inverse  is 


/ 


doa  eiujt  Jq  (b(uj2-a2)1/2)  = 


4  i  cos(a(t2-b2)1/2) 
(t2-b2)1/2 


t  <b 
t  >b 


(32) 


If  the  inverse  transform  of  3?  (  <u)/u>  can  be  found,  the  convolution 
theorem  can  be  used  once  again  to  evaluate  the  low  frequency  response 
for  eq.  (30).'  The  function  T  (t)  *  eat  S(t),  shown  in  Fig.  1,  bears  a 
reasonable  resemblance  to  the  temporal  behavior  of  an  auroral  storm 
and  gives  an  imaginary  transform  for  T  (  <jj)/ <jj  of  the  form 


Im 


/; 


d*  e1^  *  il 

U)  a 


•  at 


(33) 


The  convolution  theorem  subsequently  leads  to: 


,  '  Q  aig  e 

{> .  (x,  z,  t)  =  - 2 — 

Po(0)  a 


z/ZH 


[  /•’/ 


$  +  x 

ds 


-a(t  -  t) 


(34) 


*  .  »  g. 
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Here  and  are  given  by 


2  \  1/2  «»s 


(x  -  17) 


(38) 


and 


*—*•»)*) 


2  \  1/2  uj 


t 


(39) 


C  u)  '  (x-  r?-  cr ) 

The  remaining  integration  over  r  must  be.  carried  out  numerically. 
Velocity  perturbation: 

Because  the  ts  transform  of'the  horizontal  velocity  perturbation 
u  is  related  to  the  transform  of  the  pressure  perturbation  by  the  mo- 
mentum  equation  (5): 


U(X,  Z,U>)  =  - 


1 

its 


[  Hg 


a! 

d* 


-  F 


J 


(40) 


a  differential  equation  for  u  analogous  to  eq.  (10)  can  be  written: 


£  E  +  a  dzs 


a* 


A,  ?-=-  +  A,  E  *  M, 

a  x2  2  ^ 


2  '  '“1  n  L2 


(41) 


where 


Z  .  e-z/2H  u 


(42) 
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and 


M„ 


Hg  QUX 

im  Q  * 


for  the  domain  outside  of  the  range  77  <x  <77+  <s  directly  beneath 
the  electrojet.  Proceeding  as  in  the  calculation  of  the  pressure  per 
turbation,  a  solution  is  obtained  by  the  use  of  Fourier  transforms; 


Q'«  e*'2H 

u(x,z,t)  = - * 


/>  (0)  a* 


£  J  d  r  |  £  e"®^  "  -  2  a  (t  -  t)  £  cos  (£  +  X.  -  z) 


-cos  Bj  (|  -  z)  J  /  (x  -  77)  / <Ug^  J 
-  j"  d.r  |  £  ea(t  -  1* )  _  (t  .  tj  j  £  cos  B2(^  +  X-z) 


-  cos  B2(  £  -  *)  J  /  ^ 


B.  (x  -  t/-  a ) 


':)  I] 


where 


O  '  *  Q  C2  ./y 


The  corresponding  vertical  shear  of  the  horizontal  wind  is  also  of 
interest: 


COMPUTATIONAL  RESULTS  AND  DISCUSSION 


Figure  2  shows  some  time  histories  of  the  pressure  pertur¬ 
bation  at  various  ranges  and  altitudes  below  the  electrojet  calcul¬ 
ated  from  eq.  (37)  with  the  representative  values  A  ok  =  10^  Amps, 
a  -  12  km,  *7=0  km,  £  =  100  km,  X  =  20  km,  H  =  104  m,  g  =  9.5  m 
sec-2,  V  =  1.4,  wA=  1/2  ^Yg/H  ,  Ug  =  .017  sec-1,  B  =  6  x  10"5 
Wm  2,  and  Pq  (0)  =  J . m  for  a  =  (JO  min)-1  and  a  =  (300  min) 
Although  the  actual  behavior  and  spatial  morphology  of  the 

electrojet  are  not  we* U‘»»n  (Francis,  1974),  these  parameters  together 
with  the  model  shewn  in  Fig.  1  approximate  some  of  the  characteristics 
deduced  from  radar  observations  (Greenwald  et  al. ,  1973,  1975,  Brekke 
1978).  Since  the  mathematical  treatment  was  carried  out  for  the  limit 
of  low  frequencies, temporal  structures  with  time  scales  >10  min  are 
presumably  valid  results.  The  fast  transients  seen  in  Fig.  2  are 
questionable  features  which  may  be  related  in  part  to  high  frequency 
components  in  the  electrojet  time  structure,  especially  at  onset  (see  Fig.  1). 
Caution  must  also  be  exercised  in  the  application  of  the  quasi-linear 
calculation  at  altitudes  where  the  magnitude  of  the  perturbation  becomes 
comparable  to  the  ambient  pressure.  Because  the  solution  is  symmetric 
about  z  =  |  +  \/2,  the  oscillations  at  z  =  10  km  can  be  compared  with 
the  characteristics  of  TIDs  which  are  observed  at  an  altitude  of  200  km 
where  the  pressure  is  ~10-^  mb.  Chimonas  and  Hines  (1970)  imply  that 
typical  —  1  hr  period  perturbations  of  a  few  percent  are  seen  3000  km 
equatorward  of  the  auroral  zone,  consistent  with  the  results  shown  in 
Fig.  2. 
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PRESSURE  PERTURBATION  lantbN  ioW 


Another  view  of  the  spatial  and  temporal  development  of  the  pres¬ 
sure  perturbation  is  seen  in  Fig.  3  which  shows  how  the  sinusoidal 
oscillations,  which  belong  to  the  class  of  gravity  waves,  develop  at 
distances  from  the  source  where  significant  dispersion  of  the  Fourier 
components  occurs.  This  picture  also  illustrates  the  direction  of  phase 
propagation,  which7  is  normal  to  the  direction  of  energy  propagation 
from  the  source  as  expected  for  gravity  waves  (Beer,  1978).  These 
patterns  are  symmetric  about  z  =  £  +  \/Z  and  change  sign  at  the 
source  midplane  x  =  +  where  the  Lorentz  force  switches  direction 

between  toward  and  away  from  the  electrojet. 

An  amplified  view  of  the  very  near-f;eld  of  the  pressure  pertur¬ 
bation  is  given  in  Fig.  4.  Here,  the  maximum  disturbance,  which 
becomes  infinite  at  the  site  of  the  N-S  boundaries  of  the  model  elec¬ 
trojet  (x  -77  a  0  and  x  -  77  =  a  )  is  seen.  Because  the  sharp  vertical- 

C 

edge  geometry  of  the  model  is  unrealistic,  this  behavior  is  an  artifact; 
however,  the  nearby  perturbations  of  <  0. 1  -1 . 0  mb  may  be  realized  in 
nature  if  the  electrojet  edge  is  fairly  localized  in  range.  Of  course,  a 
perturbation  of  this  magnitude  is  too  large  to  be  considered  linearly  at 
altitudes  >  50  km. 

The  earth  reflected  wave  discussed  by  Francis  (1974)  can  be  easily 
included  if  it  is  assumed  to  be  reflected  with  efficiency  <  and  is  describ- 
able  by  an  image  source  at  z  *  -  ( £  +  \/2): 

P^tal  *  PL  (x,z,  t,  \)  +  €  px  (x,z,t,  -  X) 

An  example  of  the  modification  of  the  pressure  pattern  with  the  total 

(  «*  1)  reflected  component  added  is  shown  in  Fig.  5.  It  can  be  seen  from 
Fig.  5  that  the  pressure  pattern  becomes  more  comolicated.  but  more 
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x  =  10  km 


500  km 


1000  km 


Fig.  5.  Calculated  Pressure  Pattern  with  Earth-Reflected  Component 
Added  for  the  Case  or  =  (30  min)*l,  t  =  60  min.  Format  of 
this  contour  diagram  is  described  in  Fig.  2  caption. 


importantly,  that  the  very  near-field  magnitude  is  evidently  reduced 
substantially  by  destructive  interference  with  the  direct  perturba¬ 
tion.  Since  the  reflection  properties  of  the  boundary  layer  are  pro¬ 
bably  variable.  Fig.  3  and  Fig.  5  represent  extreme  cases. 

The  velocity  perturbation  described  by  eq.  (44)  is  quite  struc¬ 
tured  compared  to  the  pressure  perturbation  and  so  is  not  easily  dis- 
played  as  a  contour  diagram.  Figure  6  shows  some  representative 
altitude  profiles  of  the  velocity  perturbation  at  two  different  ranges. 
Examples  of  the  time  histories  of  these  perturbations  for  the  case 
a  =  (300  min)  *  are  given  in  Fig.  7.  As  noted  in  the  discussion  of 
the  pressure  perturbations  that  were  shown  in  Fig.  2,  rapid  time  scale 
(<  10  min)  fluctuations,  which  appear  to  dominate  the  velocity  spectrum 
at  x  =  20  km  where  very  large  perturbations  are  seen,  are  of  question¬ 
able  validity.  However,  low  frequency  (w<  (jg)  variations  of  ~ 20-30  m  s~* 
are  clearly  present  in  some  of  the  waveforms,  especially  at  t  «  100  s*n. 

A  more  precise  determination  of  the  low  frequency  content  ,‘:t  these 
results  could  in  principle  be  obtained  by  a  numerical  Fourier  analysis  / 
but  this  additional  computation  has  not  been  carried  out  at  the  present  time. 
The  corresponding  vertical  shears  of  the  velocity  can  be  inferred  from 
Fig.  6  or  calculated  from  eq.  (46).  The  latter  method  produced  the  altitude 
profiles  and  time  histories  displayed  in  Figs.  8  and  9,  respectively.  Even 
if  the  high  frequency  component  is  disregarded,  it  is  apparent  that  shears 
of  up  to  'vlO  m  s”* /km  occur  when  the  value  cr  =  (300  min)”*  is  used. 

Not  surprisingly,  the  long-lived  auroral  activity  produces  the  stronger 
low  frequency  effects. 
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VELOCITY  PERTURBATION  ulms*1) 


6.  Altitude  Profiles  of  Horizontal  Wind 
Velocity  Perturbation  by  Model 
Electrojet  at  Ranges  x  =  20  km  and 
x  *  100  km  Calculated  from  eq.  (44) 
for  t  =  60  min,  a  =  (30  min)“l  (a) 
and  a  =  (300  min)*1  (b).  Altitude 
resolution  is  0.  5  km.  Structure  at 
~20  km  altitude  intervals  results 
from  altitude  structure  of  source. 
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8.  Altitude  Profiles  of  Horizontal  Wind  Shear 
Corresponding  to  Cases  Shown  in  Fig.  6 
as  Calculated  from  eq.  (46): 

(a)  a  =  (30  min)-*;  (b)  a  =  (300  min)“l. 
Spatial  resolution  is  0.  5  km.  • 


SHEAR  PERTURBATION  du/dz  (variable  scale) 

106  km  z  =  78  km  z»50km 


x  =  20  km 


a  »  300  min"^ 


x=  100  km 


x  =  1000  km 


1  ms  /km 
SCALE 


-vW 

100  ms"^/km 

h*— 3W)  min— < 
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Fig.  9.  Time  Histories  of  Shear  Perturbations  at  Several 
Ranges.  Temporal  resolution  is  5  min. 
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At  this  point  in  the  discussion,  it  is  useful  to  point  out  that  all 
of  the  above  results  can  be  rather  easily  extended  to  other  electrojet 
current  strengths  and  scale  sizes.  The  former  modification  involves 
straightforward  renormalization.  Structured  or  latitudinally  wider 
electrojet  geometries  produce  perturbations  that  can  be  constructed 
from  superpositions  of  the  given  patterns,  appropriately  shifted  in 
latitude  in  accordance  with  the  source  distribution. 

Among  the  results,  velocity  shears  are  of  particular  interest 
because  they  have  been  found  to  influence  the  potential  for  air  turbulence 
at  a  specific  location.  In  particular, a  standard  index  which  is  used  as  a 
measure  of  instability  is  the  gradient  Richardson  number  (Woods,  1969) 


Turbulence  in  both  the  ocean  and  the  atmosphere  is  found  to  be  correlated 
with  large  vertical  shears  in  velocity  and  with  small  Richardson  numbers, 
but  the  definition  of  the  threshold  or  critical  values  appears  to  vary 
with  the  circumstances.  For  example,  Waco  (1970)  found  empirically 
that  stratospheric  clear  air  turbulence  (CAT)  was  likely  to  occur  in 
regions  where  the  wind  shear  exceeded  ~2.5  m  s  *  /km  and  where  Ri 
was  less  then  15.  Unfortunately,  a  similar  analysis  for  the  upper 
middle  atmosphere  does  not  appear  to  have  been  carried  out.  Neverthe¬ 
less,  it  is  interesting  to  note  that  wind  shears  with'  a  low  frequency  com¬ 
ponent  of  the  order  of  the  forementioned  magnitude  are  found  for  the  case 
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of  the  slowly  decaying  electrojet  event  if  the  current  is  reduced  from 
106  A  to  105  A.  If  one  conservatively  adopts  the  velocity  shear  profile 
at  x  =  100  km  from  Fig.  8b,  for  which  the  low  frequency  content  is 
clearly  distinguishable  in  Fig.  9,  Richardson  numbers  for  a  particular 
time  and  location  in  the  perturbed  zone  can  be  estimated.  The  average 
magnitude  of  the  Brunt- Vaisala  frequency  in  the  low  and  middle 
atmosphere  is  shown  in  Fig.  10a.  As  seen  here,  the  actual  value  of 

o>  is  not  constant,  as  assumed  in  the  calculations,  but  varies  with 
© 

altitude  due  to  the  temperature  structure  of  the  atmosphere  according 
to  the  formula 


0)  2  = 

g 


-g(x 


8  P 

8  z 


"yrt) 


(47) 


where  R  is  the  gas  constant 

If  it  is  assumed  that  the  auroral  activity  affects  the  wind  shear  more  than 

it  affects  U  ,  the  Richardson  number  in  the  vertical  cross-section 
g 

corresponding  to  the  shear  perturbation  of  Figs.  8b  and  10b  is  given  by 
Fig.  10c.  Layers  of  small  Richardson  number  occur  near  the  altitude 
of  the  electrojet,  suggesting  the  possibility  of  upper  mesosphere  turbulence 
generation  by  auroral  activity.  Although  the  effects  discussed  above  are 
restricted  to  high  altitudes,  it  is  worthwhile  to  note  that  the  local  minima 
in  < j  at  the  mesopause  and  tropopause  make  these  regions  particularly 

9 

susceptible  to  turbulence  because  Ri  for  a  particular  wind  shear  is 
automatically  smaller  there.  The  addition  of  an  ambient  wind  shear  to  the 
wave-induced  shear  could  further  increase  the  probability  of  CAT  or 
alternatively  suppress  instability  depending  on  whether  the  auroral  per¬ 
turbations  enhance  or  reduce  the  local  background  shear.  Of  course,  these 

inferences  are  based  on  a  calculation  with  many  potential  sources  of  error; 
moreover,  is  probably,  affected  by  both  background  and  aurorally  induced 
temperature  and  density  variations. 
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Fig.  10.  (a)  Approximate  Altitude  Dependence  of  Square  of  Brunt-Vaisala 

Frequency,  (b)  Velocity  Shear  Perturbation  of  Fig.  8b,  and  (c) 
Value  of  Ri  Calculated  from  these  Quantities. 


QUALIFICATIONS  AND 
CONCLUDING  REMARKS 


A  major  point  of  concern  in  the  above  calculations  is  the 
assumption  of  an  isothermal  atmosphere.  However,  Francis  (1974) 
has  argued  that  no  strong  ducting  mechanisms  operate  for  the  waves 
considered  in  both  his  own  report  and  here.  Free  propagation  below 
100  km  altitude,  with  the  exception  of  the  boundary  at  the  surface  of 
the  earth,  was  also  assumed  in  that  author's  earlier  analysis  of 
medium  scale  TIDs.  The  full  nonisothermal  calculation  must  be 
carried  out  by  methods  (e.g.  see  Friedman,  1966)  which  are  outside 
of  the  scope  of  this  work.  The  effects  of  ambient  winds  (e.g.  see  Hines 
and  Reddy,  1967)  are  also  not  incorporated  in  the  present  effort  for 
similar  reasons.  Moreover,  the  application  of  linear  perturbation 
theory  here  can  be  questioned,  as  discussed  above,  as  can  the  credibility 
of  the  high  frequency  component  of  the  calculated  results.  The  magnitude 
and  geometry  of  the  assumed  electrojet  model  may  also  be  challenged. 

The  sharp  edges  of  the  volume  containing  the  current  can  produce  un¬ 
realistic  features  in  the  solution;  furthermore,  substorms  are  often 
found  to  have  electrojet  currents  of  only  10^  A  and  meridional  extents  of 
several  hundred  km.  Thus,  the  calculated  effects  are  likely  to  over¬ 
estimate  the  usual  situation.  However,  in  spite  of  these  qualifications, 
the  present  calculations  seem  to  raise  some  interesting  issues  concerning 
the  near-field  atmospheric  perturbations  caused  by  the  auroral  electrojet 
such  as  the  possibility  of  air  turbulence  generation  in  the  auroral  zone 
during  severe  geomagnetic  storms. 


^  v 


35 


REFERENCES 


Beer,  T. ,  Atmospheric  Waves,  John  Wiley  &  Sons,  New  York,  1974. 

Blumen,  W.  and  R.  G.  Hendl,  On  the  role  of  Joule  Heating  as  a 
source  of  gravity  wave  energy  above  100  km,  J.  Atmos.  Sci.  26, 

210,  1969. 

Brekke,  A. ,  On  the  relative  importance  of  Joule  heating  and  the  Lorentz 
force  in  generating  atmospheric  gravity  waves  and  infrasound  waves 
in  the  auroral  electrojets,  J.  Atmos.  Terr.  Phys.  41^  475,  1979. 

Brekke,  A.  And  C.  L.  Rino,  High- resolution  altitude  profiles  of  the 
auroral  zone  energy  dissipation  due  to  ionospheric  currents,  J. 
Geophys.  Res.  83,  2517,  1978. 

Chimonas,  G. ,  Infrasonic  waves  generated  by  auroral  currents.  Planet. 
Space  Sci.  1_8,  591,  1970. 

Chimonas,  G.  and  C.  O.  Hines,  Atmospheric  gravity  waves  launched 
by  auroral  currents,  Planet.  Space  Sci.  Jl_8,  565,  1970. 

Chimonas,  G.  and  W.  R.  Peltier,  The  bow  wave  generated  by  an  auroral 
arc  in  supersonic  motion.  Planet.  Space  Sci.  JJ5,  599,  1970. 

Chiu,  Y.  T. ,  Planetary  scale  wave  response  to  auroral  heating  of  the 
neutral  atmosphere,  J.  Geophys.  Res.  _81,  1231,  1976. 

Davis,  M.  J. ,  On  polar  substorms  as  the  source  of  large-scale  traveling 
ionospheric  disturbances,  J.  Geophys.  Res.  76,  4525,  1971. 

37  —  , 


FK&CLDl.tQ  tAQ K  E 


Francis,  S.  H.,  A  theory  of  medium-scale  traveling  ionospheric 
disturbances,  J.  Geophys.  Res.  79,  5245,  1974. 

Friedman,  J.  P. ,  Propagation  of  internal  gravity  waves  in  a 

thermally  stratified  atmosphere,  J.  Geophys.  Res.  7^,  1033,  1966. 

» 

Greenwald,  R.  A.,  W.  L.  Ecklund  and  B.  B.  Balsley,  Auroral 
currents,  irregularities  and  luminosity,  J.  Geophys.  Res.  78_, 

8193,  1973. 

Hines,  C.  O.  and  C.  A.  Reddy,  On  the  propagation  of  atmospheric 
gravity  waves  through  regions  of  wind  shear,  J.  Geophys.  Res. 

72,  1015,  1967. 

Hunsucker,  R.  D. ,  Estimate  of  the  relative  importance  of  Joule  heating 
and  the  Lorentz  force  in  generating  atmospheric  gravity  waves  from 
the  auroral  electrojet,  J.  Geophys.  Res.  82,  4826,  1977. 

Richmond,  A.  O.  and  S.  Matsushita,  Thermospheric  response  to  a 
Magnetic  Substorm",  J.  Geophys.  Res.  80,  2839,  1975. 

Testud,  J. ,  Gravity  waves  generated  during  magnetic  substorms, 

J.  Atmos.  Terr.  Phys.  32,  1793,  1970. 

Waco,  D.  E.,  A  statistical  analysis  of  wind  and  temperature  variables 
associated  with  high  altitude  clear  air  turbulence  (HICAT),  J.  Appl. 
Meteor.  9,  300,  1970. 


Wilson,  C.  R.,  Infrasonic  wave  generation  by  aurora,  J.  Atmos. 
Terr.  Phys.  37,  973,  1975. 

Woods,  J.  D.,  On  Richardson's  number  as  a  criterion  for  laminar- 
turbulent -laminar  transition  in  the  ocean  and  atmosphere.  Radio 


LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  la  conducting  exper¬ 
imental  and  theoretical  Investigations  necessary  for  the  evaluation  and  applies' 
tlon  of  scientific  advances  to  new  military  space  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  In 
dealing  with  the  many  problems  encountered  in  the  nation's  rapidly  developing 
space  systems.  Expertise  in  the  latest  scientific  developments  Is  vital  to  the 
accomplishment  of  tasks  related  to  these  problems.  The  laboratories  that  con' 
tribute  to  this  research  are: 

Aerophy8lcs  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer,  propulsion  chemistry  and  fluid  mechanics,  structural  mechanics,  flight 
dynamics;  high'teaperature  thermomechanics,  gas  kinetics  and  radiation;  research 
in  environmental  chemistry  and  contamination;  cv  and  pulsed  chemical  laser 
development  including  chemical  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo- 
spheric  optics^  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  In  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenomena,  thermionic  emission,  photosensitive  materials 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmental  research  and 
monitoring. 

Electronics  Research  Laboratory:  Microelectronics,  CaAs  low-noise  and 
power  devices,  semiconductor  latere,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  comauni cat ions,  Ildar,  and  electro-optics; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  imaging;  millimeter-wave  and  microwave  technology. 

Information  Sciences  Research  Office:  Program  verification,  program  trans¬ 
lation  per  f  ormance- sens 1 1 lve  system  design,  distributed  architectures  for 
spaceborne  computers,  fault-tolerant  computer  systems,  artificial  intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  Forms  of  carbon;  component  failure  analysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  In 
apace  environment;  materials  performance  In  space  transportation  systems;  anal¬ 
ysis  of  systems  vulnerability  and  survivability  In  enemy-induced  environments. 

Space  Sciences  Laboratory:  Atmospheric  and  Ionospheric  physics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglow;  magnetospherlc  physics,  cosmic  rays,  generation  and  propagation  of 
plasma  waves  in  the  magnetosphere;  solar  physics,  infrared  astronomy;  the 
affects  of  nuclear  axploslons,  magnetic  storms,  and  solar  activity  on  the 
earth's  atmosphere,  Ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  end  particulate  radiations  in  space  on  apace  systems. 


